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The Inert Doublet Model is an extension of the Standard Model including one extra "Inert 
scalar doublet" and an exact Z2 symmetry. The "Inert scalar" provides a new candidate for 
dark matter. We present a systematic analysis of the dark matter abundance assuming the 
standard freeze-out mechanism and investigate the potentialities for direct and gamma indirect 
detection. We show that the dark matter candidate saturates the WMAP dark matter density 
in two rather separate mass ranges, one between 40 and 80 GeV, the other one over 400 GeV. 
We also show that the model should be within the range of future experiments, like GLAST 
and EDELWEISS II or ZEPLIN. 



1 Introduction 



Many evidences for the existence in the universe of dark matter has been put forward over the 
years. It can be inferred from the dynamics of galaxies and of clusters of galaxies, from analysis 
of the CMBR, from structure formation, etc. The question that arises then is what is the nature 
of dark matter and which extension of the Standard Model do we have to consider in order to 
account for these observations? A profusion of models of dark particles have been proposed over 
the years and it is much hoped that present and forthcoming experiments will throw some light 
on the matter (for a review, see for instance^). 

In these proceedings, we study a simple extension of the Standard Model with one extra 
scalar doublet and an exact Z2 symmetry. In this framework, the candidate for dark matter 
is one of the two neutral scalars arising from the extra doublet. The latter was called "Inert 
doublet" by Barbieri et al in^ because it has no direct coupling to matter fields. However it 
couples to the standard gauge fields. The phenomenology of its neutral and charged components 
is quite simple and yet very rich. 

This is not the only attractive feature of the model. As was pointed out in™, the Inert Doublet 
Model (IDM) could allow for a Higgs mass up to 500 GeV still fulfilling the LEP Electroweak 



Precision Test measurements. Here, we will only consider a Higgs with a mass of 120 GeV. In the 
reference paper^, one can find a more detailed study of the IDM with Higgs masses of 120 GeV 
and 200 GeV (See also - where the authors considered Higgs masses up to 500 GeV). Another 
interesting aspect of models like the IDM is that it could pave the way toward an understanding 
of the relation between the abundance of dark and ordinary matter (see e.g. ^). 



2 Short description of the Model 

The IDM is a particular two Higgs doublet model, in which one of the doublet, Hi, plays the 
role of the standard Brout-Englert-Higgs doublet while the second one, H 2 , is the source for dark 
matter candidates. In order to guarantee the stability of the dark matter particles, one invoke a 
Z 2 symmetry under which all Standard Model fields are even and 

H x -> H x and H 2 -» -H 2 . 

This discrete symmetry also prevents the appearance of flavor changing neutral currents in 
this model. Moreover, we assume that Z 2 is not spontaneously broken. This model was first 
introduced by Deshpande and 

Mal%seealsc£3), 

and the dark matter aspect was recently discussed 
by Cirelli et afe and Barbieri et a$. Their initial purpose and some of their assumptions were 
nevertheless not exactly identical. In addition, the neutral scalar reaching the dark matter 
WMAP abundance was found to be in the mass range of 60 to 75 GeV for Barbieri et while 
for Cirelli et av^ it was of order of 430 GeV. We first study the details of the model before to 
elucidate this apparent incompatibility. 

The most general potential of the model can be written as 

V = fil\H 1 \ 2 + f il\H 2 \ 2 +X 1 \H 1 \ 4 +X 2 \H 2 \ 4 +X 3 \H 1 \ 2 \H 2 \ 2 +X A \HlH 2 \ 2 + ^ \(h\h 2 ) 2 + h.c. 



The vacuum expectation value of Hi is given by (Hi) = -^g with v = \J — [i 2 /X\ = 248 GeV, 

while assuming for simplicity fi 2 > 0, we have (H 2 ) = 0. The mass of the Higgs particle h 
is M 2 = — = 2Xiv 2 while the mass of the charged, H + , and two neutral, #0 and Aq, 
components of the field H 2 are given by 



M 2 H+ = f4 + X 3 v 2 /2 



M 2 Ho = f4 + (A3 + A 4 + X 5 )v 2 /2 

M 2 Ao = (4 + (X 3 + X 4 -X 5 )v 2 /2. (2) 



For appropriate quartic couplings, Hq or Aq is the lightest component of the H 2 doublet. 

In the absence of any other lighter Z2-odd field, either one is a candidate for dark matter. For 

definiteness we choose Hn. All our conclusions are unchanged if the dark matter candidate is An 
[oi 

instead. Following 1 ^ 1 , we parameterize the contribution from symmetry breaking to the mass of 
Ho by Al = (A3 + A4 + \§)/2, which is also the coupling constant between the Higgs field h and 
our dark matter candidate Hq. 



3 Dark matter abundance 

As inland in^, we consider a thermal production of the cold relic Hq. We have computed the 
relic abundance of Hq using micrOMEGAs2.0, a new and versatile package for the numerical 
calculation of dark matter abundance from thermal freeze-out^. 

We first present the results for fixed Inert scalars mass differences in the contour plots of 
figure [H We work in the (Mh ,^2) plane, as a result the diagonal line corresponds to Xl = 0, 
i.e. to no coupling between Hq and the Higgs boson. Away from this line, A^ increases, with 
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Figure 1: Relic density contours in the (Mh () ,H2) plane for Mu = 120. Left: low mass regime with Inert scalars 
mass differences AMAq = 10 GeV and AMH + = 50 GeV . Right: high mass regime with Inert scalars mass 

differences AMA = 5 GeV and AMH+ = 10 GeV. 



Xl < (resp. A^ > 0) above (resp. below) the diagonal. Also, we write AMAq = Ma — Mh 
and AMH C = M H + - M Ho . 

The shaded areas in the plots correspond to regions that are excluded by several constraints. 
In order not to conflict with LEP data, the mass of the H + should be larger than 79.3 GeV and 
M H o + M A o < Mz- These constraints translate into the excluded region 1 on the plot on the 
left of figure [Q The vacuum stability constraint contributes largely to the exclusion of Xl < 
couplings. This corresponds to shaded area in the domains \i2 > Mh in the plots of figure [TJ 
The remaining shaded regions are excluded due to large couplings |Aj| > in. Moreover regions 
where the couplings range as 1 < |Aj| < 47r, which is still tolerable, are shown with horizontal 
lines. The areas between two dark lines correspond to regions of the parameter space such that 
0.094 < n DM h 2 < 0.129, the range of dark matter energy densities consistent with WMAP data. 

We immediately see that there are two qualitatively distinct regimes, depending on whether 
the Hq is lighter than the W and Z and/or the Higgs boson. For the low mass regime, let us study 
figureCH the plot on the left. The two processes relevant below the W, Z or h threshold are the Hq 
annihilation through the Higgs and Hq coannihilation with Aq through Z exchang^]. Both give 
fermion-antifermion pairs, the former predominantly into bb. Coannihilation into a Z may occur 
provided AMAq is not too important, roughly AMAq must be of order of Tj Q ~ Mh /25. As 
the mass of Hq goes above W, Z or h threshold, Hq annihilation into WW, ZZ and hh become 
increasingly efficient, an effect which strongly suppresses the Hq relic density. The region 4 of 
figure Q] , corresponding to Mh G [40, 80] GeV, appears to be the only region consistent with 
WMAP data. 

For the high mass regime, we can derive the general trends from figure [TJ the plot on the right. 
No new annihilation channel opens if Mh is heavier than the Higgs or the gauge bosons. There 
are then essentially two kinds of processes which control both the abundance: the annihilation 
into two gauge bosons, dominant if ^2 < Mh , and the annihilation into two Higgs, which 
dominates if [12 > Mh . Coannihilation plays little role. 

The abundance of dark matter is suppressed over most of the area of the plot because of 
large quartic coupling effects on the cross-sections. Let us emphasize that in this regime, it is 



a HoH + coannihilation is suppressed for our choice of AMH C . 




Figure 2: Left: Relic density and Right: scattering cross-section intervening in direct detection searches, all as 
a function of the mass of dark matter and comparison with the MSSM. For the direct detection plot, the light 
colors correspond to 0.01 < QdmIi 2 < 0.3, while the dark colors correspond to 0.094 < SIdm/i 2 < 0.129. 



possible to reach agreement with WMAP data, but only at the price of some fine tuning. We 
need to keep the mass splittings between the components of H2 relatively small. First because 
large mass splittings correspond to large couplings and second because the different contributions 
to the annihilation cross-section must be suppressed at the same location, around Al = (i.e. 
Mh — A 1 — Ma ~ Mh + in this case). As it can be seen in figure [U the plot on the right, 
the area consistent with WMAP corresponds to the narrow region around the diagonal with 
M H > 800 GeV for AMA = 5 GeV and AMH C = 10 GeV. Notice that this behavior is limited 
by the unitarity bound on the total annihilation cross-section^ which constrains the mass of 
the dark particle to be M Ho < 120 TeV. 

In figure El the plot on the left, we show a scatter plot of Qr>Mh 2 as a function of the mass 
of the dark matter candidate Mdm for a fair sample of IDMs (scanning on several Inert scalar 
mass splittings) and, for the sake of comparison, for the MSSM. We clearly see the two regimes 
(low mass and high mass) of the IDM that may give rise to a relevant relic density (i.e near 
WMAP). The MSSM models have a more continuous behavior, with O(100 GeV) dark matter 
masses. As a conclusion, the IDM provides dark matter candidates with masses as small as 40 
GeV and as large as 600 GeV in contrast with the MSSM more concentrated around ~ 100 GeV. 



4 Direct detection 



Direct detection searches look for signals of dark matter in low background detectors trying to 
measure the energy deposited by the scattering of a dark matter particle with a nucleus of the 
detector. Assuming that the main interaction contributing to the .ffo-quarks interaction is the 
spin independent H$q H$q interaction!] it can be shown^l that the Hq elastic scattering cross 
section off a proton scales like 

a Ho „ p K\ L 2 /(M Ho Ml) 2 . (3) 
In figure[21 the plot on the right, we show a scatter plot of log 10 <jdm-p as a function of Mom 

'The experiments have reached such a level of sensitivity that the Z exchange contribution Hoq Aoq is 
excluded by the current experimental limits^. Consequently, to forbid Z exchange by kinematics, the mass of the 
Ao particle must be larger than the mass of Hq by a few 100 keV. 




Figure 3: Integrated gamma-ray flux from the Galactic Center (GC) resulting from dark matter annihilation as 
a function of the mass of the dark matter candidate for the same sample of models than for direct detection 
in Fig[2] Again, the light colors correspond to 0.01 < Q.DAih 2 < 0.3, while the dark colors correspond to 
0.094 < floAih 2 < 0.129. For the plot on the left, we the took an isothermal profile (flat profile), for the plot on 
the right, we took a NFW profile (more steeper at Galactic Center). 



for the IDMs considered in the abundance plot of figure [2] that account for 0.01 < Q,DAih 2 < 0.3. 
We see that the low mass regime candidates could be detected by future experiments such as 
EDELWEISS II or by the ton sized experiments such as ZEPLIN. For the higher mass regime 
however, there is no hope for future detection in low background detector. Indeed, the WMAP 
requirement for dark matter relic density constrains the Xl couplings to be vanishing while the 
same couplings drive the amplitude of the matter- Hq scattering cross-section. 



5 Indirect detection 

The measurement of secondary particles coming from dark matter annihilation in the halo of the 
Galaxy is another promising way of deciphering the nature of dark matter. Let us emphasize that 
this possibility depends however not only on the properties of the dark matter particle, through 
its annihilation cross-sections, but also on the astrophysical assumptions made concerning the 
distribution of dark matter in the halo that supposedly surrounds our Galaxy. 

In Figure El we show the log of the produced gamma-ray flux from dark matter annihilation 
at the Galactic Center <3? 7 as a function of Mqm for the same sample of models than for direct 
detection. We computed the gamma-ray flux for the plot on the left assuming an isothermal dark 
matter density profile while for the plot on the right we assumed a Navarro, Frank and White 
(NFW) profile. The main difference between these two profiles is the slope of the dark matter 
density as a function of the galactic radius in the central part of the Galaxy. The isothermal 
profile is flat while the NFW profile is more cuspy (i.e. steeper). We see that for steeper profile 
the gamma-ray flux is larger. 

The particle physics dependence of $ 7 also clearly show up in figure 03 Indeed we see that 
( I ) 7 behaves differently in the low and the high mass regime of the IDM given that the processes 
contributing to the annihilation cross-section are different. Moreover, notice that the IDM dark 
matter candidates have typically higher detection rates than the neutralino in SUSY models, 
especially at high mass. Let us stress that the figures for indirect detection were obtained taking 
into account annihilation processes at three level only (see^, for a recent study of the IDM 



including processes at one-loop). It can be inferred from figure [3] that the IDM can give the 
right relic abundance in a range of parameters which will be probed by GLAST for NFW dark 
matter profiles. GLAST will however have no chance to observe the gamma-ray flux produced 
by annihilating Hq at the Galactic Center for flatter profiles such as the isothermal one. 

6 Conclusion 

We carried out a rather detailed analysis of the IDM as a dark matter model assuming the 
standard freeze-out mechanism. We recovered the results of Barbieri et al. and Cirelli et al. 
which a priori did not seem to match. This is because the IDM provides dark matter candidates 
in two rather separate mass ranges, one between 40 and 80 GeV, the other one over 400 GeV. 
The physics driving the existence of dark matter in these regions of the parameter space is quite 
different. 

We have also investigated the prospects for direct and indirect detection searches. Concerning 
direct detection searches, the low mass regime candidates should be detected with the futures 
ton sized experiments while the high mass regime will stay out of reach. For indirect detection 
searches we looked at the gamma-ray flux generated at the Galactic Center by dark matter 
annihilation. Whatever the dark matter density profile assumed, we have come to the conclusion 
that the Inert scalars have typically higher detection rates than the neutralino in SUSY models, 
especially at high mass. Moreover, the IDM could be probed by the future GLAST experiment. 
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